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ABSTRACT 

We describe a new method for accurately determining total gas-phase abundances 
for the Galactic halo interstellar medium with minimal ionization uncertainties. For 
sight lines toward globular clusters containing both ultraviolet-bright stars and radio 
pulsars, it is possible to measure column densities of H I and several ionization states 
of selected metals using ultraviolet absorption line measurements and of H II using 
radio dispersion measurements. By measuring the ionized hydrogen column, we min- 
imize ionization uncertainties that plague abundance measurements of Galactic halo 
gas. We apply this method for the first time to the sight line toward the globular 
cluster Messier 3 [(/, b) = (42?2, +78?7); d = 10.2 kpc, z = 10.0 kpc] using Far Ul- 
traviolet Spectroscopic Explorer and Hubble Space Telescope ultraviolet spectroscopy 
of the post-asymptotic giant branch star von Zeipel 1128 and radio observations by 
Ransom et al. of recently-discovered millisecond pulsars. The fraction of hydro- 
gen associated with ionized gas along this sight line is (45 ± 5)%, with the warm 
(T ~ 10 4 K) and hot (T > 10 5 K) ionized phases present in roughly a 5:1 ratio. 
This is the highest measured fraction of ionized hydrogen along a high-latitude pulsar 
sight line. We derive total gas-phase abundances log iV(S)/iV(H) = —4.87 ± 0.03 
and log iV(Fe)/iV(H) = —5.27 ± 0.05. Our derived sulfur abundance is in excellent 
agreement with recent solar system determinations of Asplund, Grevesse, & Sauval. 
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However, it is —0.14 dex below the solar system abundance typically adopted in stud- 
ies of the interstellar medium. The iron abundance is ~ —0.7 dex below the solar 
system abundance, consistent with the significant incorporation of iron into interstel- 
lar grains. Abundance estimates derived by simply comparing S II and Fe II to H I 
are +0.17 and +0.11 dex higher, respectively, than the abundance estimates derived 
from our refined approach. Ionization corrections to the gas-phase abundances mea- 
sured in the standard way are, therefore, significant compared with the measurement 
uncertainties along this sight line. The systematic uncertainties associated with the 
uncertain contribution to the electron column density from ionized helium could raise 
these abundances by < +0.03 dex (+7%). Uncertainties in the amount of very hot gas 
(T ~ 10 6 K) along the line of sight could also affect these determinations. 

Subject headings: ISM: atoms - ISM: structure - ultraviolet: ISM 



1. Introduction 

Studies of gas-phase elemental abundances of the diffuse interstellar medium (ISM) of the 
Milky Way and other galaxies can give important information about dust composition, interstellar 
mixing, and chemical evolution. Most absorption line studies of gas associated with the warm 
neutral medium (WNM) of the Milky Way and other galaxies derive gas-phase abundances by 
comparing columns of the dominant ionization stage of a metal species in the WNM 1 with the 
column of neutral hydrogen, assuming (often implicitly) no ionization corrections are needed. 
Sight lines through the multiphase ISM of a galaxy not only probe neutral material, but often 
include gas associated with the warm ionized medium (WIM) and hot ionized medium (HIM). 
Ionized gas can represent up to ~ 50% of the total hydrogen column density for sight lines through 
the Galactic halo, with an average value of ~ 25% (Reynolds 1993). While the WIM likely 
contains very little H I, the dominant ionization states of many metals are the same for gas in the 
WNM and the WIM (e.g., S II, Fe II, Si II - see Sembach et al. 2000 and Haffner, Reynolds, & 
Tufte 1999). Thus, comparing the column densities of the dominant ionization states of metals to 
the column density of neutral hydrogen along a sight line will give apparent metal abundances that 
are too high compared with the true values. The neglect of ionization corrections may be a source 
of non-negligible systematic uncertainties to measurements of interstellar gas-phase abundances 
(Sembach et al. 2000; Howk & Sembach 1999). 



We refer to the dominant ionization states of elements in the warm neutral medium as "low ions." The dominant 
ionization stage of an element is typically the first whose ionization potential is higher than that of neutral hydrogen. 
We refer to the next higher ionization states as "intermediate ions." 
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In this paper we present a new approach for studying gas-phase abundances along extended 
paths through the multiphase halo of the Milky Way that circumvents such ionization uncertainties. 
Our strategy is to observe ultraviolet (UV) absorption lines from low and intermediate ions toward 
UV-bright stars in globular clusters that contain pulsars. This allows us to measure all of the 
ionization states of certain metals (e.g., sulfur and iron) present in the WNM and WIM in these 
directions. The total hydrogen column to which these metals are compared is the sum of the H I 
and H 2 column densities, derived from UV spectra, and the H II column density, determined using 
radio observations of the dispersion measures to the pulsars (which gives the electron column 
density). We present a full discussion of our methodology in §2. 

We use this method to study the gas-phase abundances along the sight line toward the post- 
asymptotic giant branch star von Zeipel 1128 (vZ 1128) in the distant globular cluster Messier 
3 (M 3). This sight line probes gas in the first 10 kpc above the Galactic plane. We discuss 
the UV absorption line and pulsar dispersion measure observations of this sight line in §3. The 
general properties of this interstellar sight line have previously been discussed by Howk, Sembach, 
& Savage (2003; hereafter Paper I). We also apply our new strategy for determining gas-phase 
abundances in this section, showing that the neglect of ionization corrections can lead to significant 
errors in the derived metal abundances along this (and presumably other) high-latitude sight lines. 
We end with a general discussion of the implications of this work in §4. 



2. Methodology 

In a multiphase medium, the total abundance, A(X), of an element X with respect to hydro- 
gen is 

A(X) = iV(H I) + iV(H II) + 2iV(H 2 ) ' (1) 

where N(X j ) is the column density of the jth ionization stage of X. The numerator represents 
the total column density of the metal X, while the denominator is iV(H I) + N(H II) + 2iV(H 2 ) = 
iV(H), the total hydrogen column density. It is unusual that all the terms in both the numerator 
and the denominator can be measured. One approach is to identify special sight lines along which 
this is possible, such as those toward globular clusters that contain both UV-bright stars and radio 
pulsars. 

The derivation of metal ion column densities is done using UV absorption line measurements 
against background sources (in our case against UV-bright stars) using the Far Ultraviolet Spec- 
troscopic Explorer (FUSE) and/or the UV spectrographs on-board the Hubble Space Telescope 
(HST). A review of UV absorption line abundance techniques and measurements is given by Sav- 
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age & Sembach (1996). The transitions accessible to these instruments (with wavelength ranges 
of 912 < A < 1187 A for FUSE and 1150 < A < 3100 A for HST) limit the species for which 
we can reasonably approximate the total summation in the numerator of equation (1). Sulfur and 
iron are the only elements for which measurements are typically feasible for all of the ionization 
states expected to contribute significantly in the WNM and WIM (see Table 1 of Howk & Savage 
1999 for a summary of WIM tracers). Measurements of gas associated with the HIM can also be 
made using absorption line spectroscopy of the O VI 1031.926 and 1037.617 A doublet. However, 
it is not possible to measure the ions of S and Fe that dominate in the HIM. For example, the dom- 
inant ionization states of sulfur for logT > 5.5 are S VII and higher (Sutherland & Dopita 1993), 
which have no measurable UV transitions. We will make a correction (described below) for these 
still-unseen ions when deriving A(X). 

The neutral and molecular hydrogen column densities used in the denominator of equation 
(1) are also measured using UV absorption line techniques and instruments. H I column densities 
along sight lines to UV background sources are usually measured using the strong damping wings 
of the hydrogen Lyman-a transition at 1215.670 A (e.g., lenkins 1971; Bohlin, Savage, & Drake 
1978; Diplas & Savage 1994). Measurements of H 2 are made using absorption lines of the Lyman 
and Werner electronic transitions in the FUSE bandpass (e.g., Savage et al. 1977; Shull et al. 
2000). 

We derive the H II column density from the radio dispersion measure, DM = N(e~) = 
J n e dl, observed toward pulsars in the target globular clusters. We make two corrections to the 
observed N(e~) in order to derive iV(H II) from N(e~): (1) a correction for the contribution to 
the electrons from ionized helium (the contribution from heavier elements is negligible); and (2) 
the removal of the contribution of the HIM, for which metal ions are typically not measurable. 

The electron column in H+He gas is N(e~) = iV(H II) + iV(He II) + 2iV(He III). We define 
the helium correction factor rj such that iV(H II) = r]N(e~). Therefore, 



1 + A(He)a;(He + )/:r(H+) + 2A(He)a;(He +2 )/x(H+) ' 

where the ionization fraction of an ion X J in the ionized gas is x(X^) = N(Xi)/N(X). We 
assume A(Re) = 0.1 throughout. Values for rj will typically be in the range 0.8 to 1.0. 

Emission line studies of the WIM indicate 0.67 < x(R + ) WIM < 1.0 (Reynolds 1989; 
Reynolds et al. 1998) and x(Re + ) WIM < x(H. + ) WIM (Reynolds & Tufte 1995). Most pho- 
toionization models of the WIM predict x(Re +2 ) WIM to be very small (Sembach et al. 2000; 
Mathis 2000), giving values of r] WIM near 1.0. However, Arabadjis & Bregman (1999) suggest 
x{Rq +2 )wim is non-negligible on the basis of X-ray absorption estimates (although see Slavin, 
McKee, & Hollenbach 2000 for alternative interpretations of these data). In what follows, we con- 
sider two cases: x(Re +2 ) WIM ~ and x(Yie +2 ) WIM = 1 — x(Re + ) WIM ', we refer to these as the 
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minimum and maximum helium ionization cases, respectively, and tabulate characteristic values 
of 77 in Table 1 . We favor the minimum helium ionization model given that much of the sulfur in 
the WIM is in the form of S II (Haffner et al. 1999), implying the ionization state of metals in the 
WIM is relatively low. 

It is extremely difficult to measure the dominant metal ionization states in the HIM. We there- 
fore subtract an estimated HIM contribution, N(e~) HIM , from N(e~) to derive a value of iV(H II) 
appropriate for the WIM alone. We use measurements of O VI to estimate N{e~)niM, adopting 
the following two-component prescription: 



iV(e-) 



HIM — 



Vhim 



iV(OVl) . . ' 

+ (n H ) c d 



A(0)x(0+ 5 ) 



(3) 



where iV(0 VI) is the measured O VI column density, and d is the distance of the globular cluster 
from the sun. The first term in brackets is the column density of T ~ 3 x 10 5 K gas associated with 
the thick disk distribution of the Milky Way, which is best probed by O VI (Wakker et al. 2003; 
Savage et al. 2003). The second term accounts for an extended hot corona assumed to be uniformly 
distributed and characterized by T > 10 6 K and a total hydrogen density (nn) c ~2x 10~ 4 cm -3 . 
The existence of this component is derived from observations of O VI associated with high- velocity 
clouds, which suggest these clouds are interacting with a hot, low-density plasma extending to large 
distances (> 70 kpc) from the Sun (Sembach et al. 2003). Supporting evidence for a hot corona 
is provided by very strong zero redshift X-ray absorption lines of Ne IX, O VII and O VIII in the 
spectra of X-ray bright active galaxies (Nicastro et al. 2002; Fang et al. 2003; Rasmussen, Kahn, 
& Paerels 2003) and Galactic X-ray sources (Futamoto et al. 2004; Yao & Wang 2005). 

Defining the WIM electron column as N(e~)wiM = N(e~) — N(e~) HIM , i.e., correcting for 
the contribution from the HIM, we rewrite Eqn. (1) as 

A(x) = j(WNM,WIM) 

1 ) [N(Hl)+ VwiM N(e-) WIM + 2N(R 2 )Y 

noting the summation is over metal ionization states found in the WNM and WIM (i.e., excluding 
the HIM). We suggest using numerical values for the minimum and maximum helium ionization 
cases of rj WIM = 0.98 ± 0.01 and 0.81 ± 0.04, where the uncertainties correspond to the dispersion 
of values found in Table 1 . 

The differences in the helium ionization assumptions represent the dominant source of sys- 
tematic uncertainty in our method. The magnitude of the uncertainty is roughly 0.2 x iV(H \\) WIM /N{R), 
or < 10% for high-latitude sight lines. While spatial variations of N(e~) between the sight lines to 
the M 3 pulsars and the UV bright stars could be another potential source of systematic uncertainty, 
it seems in most cases this effect is negligible. For example, Freire et al. (2001) find the dispersion 
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in the DM values toward 15 pulsars in 47 Tuc is only 0.5% of the average with a full range of 
« 2% of the average (for pulsars with a maximum angular separation of 2'). 



3. Application to the Sight Line Toward von Zeipel 1128 in Messier 3 

We now apply the approach outlined above to determine average gas-phase abundances of 
sulfur and iron in the WNM and WIM toward the high-latitude globular cluster M 3 [(1,6) = 
(42?2, +78?7); d = 10.2 kpc, z = 10.0 kpc]. We describe in the following subsections the UV 
data and analysis that are used to derive the metal and neutral hydrogen column densities, the radio 
observations used to derive the electron column density toward M 3, and the resulting analysis of 
the total abundances of S and Fe in this direction. 



3.1. UV Absorption Line Measurements of Metals and Neutral Hydrogen 

In Paper I we presented UV absorption line measurements of the ISM in the direction of 
vZ 1128. That work was based upon FUSE observations covering the spectral range 905-1187 
A with a resolution of ~ 20 km s _1 . In this work we present new observations of vZ 1128 ob- 
tained with the Space Telescope Imaging Spectrograph (STIS) on board HST under guest observer 
proposal GO 9150. Proffitt et al. (2001) discuss the detailed characteristics of this instrument. 

The properties of our STIS observations of vZ 1128 are summarized in Table 2. A total 
of 12.7 ksec of exposure time was collected with each of the E140M and E230M intermediate 
resolution echelle gratings. These data have signal-to-noise ratios varying from ~ 10 to 20 per 
resolution element. The shorter wavelength E140M grating observations have a resolution of ~ 6.5 
km s _1 (FWHM), while the longer wavelength E230M observations have a resolution ~ 10 km s _1 
(FWHM). All observations were made through the 0'.'2 x 0'.'06 apertures. The STIS data were 
calibrated using v2.17 of the CALSTIS pipeline. 

We use the STIS data to derive column densities and limits for the species H I, S I, S II, and 
Fe I. The details of our measurements of metal line and H I column densities are presented below. 
We do not make use of the longest-wavelength E230M exposure (archive ID O6F501030). An 
analysis of the full STIS spectrum will be presented in a future work. 
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3.1.1. H I Column Density from Lyman-a 

We derive the interstellar H I column density of the vZ 1 128 sight line by fitting the damping 
wings of the Lyman-a profile. We generally follow the procedures described by Sonneborn et al. 
(2002) for determining the H I column density and uncertainty. 

Figure 1 shows the STIS spectrum of the Lyman-a absorption line toward vZ 1128. This 
absorption line contains contributions from both the stellar atmosphere of vZ 1 128 and the ISM in 
this direction. We use a model stellar atmosphere to normalize the data during the fitting process. 
Our adopted atmosphere, kindly provided by P. Chayer (2002, private communication), was cal- 
culated using the SYNSPEC package of I. Hubeny (2000, private communication) and assumed 
T eff = 35, 000 K, log g = 4.0, and He/H =0.1 (Dixon, Davidsen, & Ferguson 1994). The adopted 
atmosphere is shown as the blue line in Figure 1 . We adopt a second-order Legendre polynomial 
correction to the model atmosphere to match it to the STIS spectrum, following Sonneborn et al. 
(2002). This is done to account for any calibration and order combination uncertainties in the 
data or any uncertainties in the large-scale flux distribution of the models. The parameters of the 
polynomial are allowed to vary during the fitting process. 

The interstellar column density that best fits the STIS observations of Lyman-a is log N(H I) = 
19.98 ± 0.03. This is slightly different (0.01 dex higher) than the value quoted in Paper I and is due 
to a better calibration and combination of the data. The best-fit profile is shown as the red line in 
Figure 1 . Our error estimate includes contributions from the uncertainties in the adopted properties 
of the stellar model (c.f., Sonneborn et al. 2002), which dominate the total error budget. 

Two observations of the H I 21 cm emission line in this direction give slightly higher column 
density estimates (as discussed in Paper I). Danly et al. (1992) collected H I 21 cm emission line 
observations toward vZ 1 128 using the NRAO Green Bank 140-ft. telescope, which has a beam of 
21'. The Leiden-Dwingeloo Survey (LDS; Hartmann & Burton 1997), which has a beam size of 
35', includes a pointing centered ps 12' from vZ 1128. Table 3 gives the H I column densities in 
the direction of M 3/vZ 1128 derived from the two 21 -cm observations and the Lyman-a profile. 
The NRAO and first LDS value were derived here through a direct integration of the brightness 
temperature distribution over the range —75 < wlsr < +75 km s _1 in the optically-thin limit 
(optical depth effects are < 1% for this high-latitude sight line). The second LDS H I value is 
from Wakker et al. (2003); it is the sum of three Gaussians fit to a weighted average of four 
LDS pointings in this direction. The difference in column densities between these two treatments 
of the LDS data is likely due to the slightly different spectra adopted (a single pointing versus a 
weighted average of four pointings) and the different techniques used to derive column densities 
(a direct integration versus Gaussian fits). Both the Danly et al. (1992) and LDS observations have 
been corrected for stray radiation. The uncertainties given in Table 3 do not include contributions 
from errors in this correction. Lockman & Savage (1995) estimate the stray radiation correction 
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uncertainty of observations similar to those of Danly et al. (1992) to be of order <r[iV(H I)] ~ 10 19 
cm -2 , or ~ 0.04 dex for the column densities in Table 3. 

All of the H I columns derived using 21 -cm observations are higher than that derived above 
via a direct fit to Lyman-a absorption toward vZ 1 128. The largest difference is 0.08 dex, or 20%, 
between the NRAO spectrum and the fit to Lyman-a. This difference is likely due to structure in the 
H I column on scales smaller than the 21 cm beam (see Savage et al. 2000), although it could also 
be caused by systematic effects (e.g., unaccounted for uncertainties in the adopted stellar Lyman-a 
profile or uncertainties in the 21 cm stray radiation correction). In Figure 1 we show the best-fit 
Lyman-a profile if we assume logiV(H I) = 20.06 as the green line. The polynomial continuum 
correction is allowed to vary when fitting this profile to the data. This column over-predicts the 
optical depth in the wings of Lyman-a. In what follows, we assume the Lyman-a profile provides 
the best estimate of iV(H I) in this direction. 



3.1.2. Metal Line Column Densities 

We have used the new STIS data to derive column densities and limits for the metal species 
S I, S n, and Fe I. Our measurements of the metal line properties toward vZ 1 128 follow Paper I. 
We have fit the stellar continuum in the regions surrounding metal transitions using low order 
Legendre polynomials. Following Sembach & Savage (1992), we have directly integrated the 
observed line profiles and apparent optical depths to arrive at the equivalent widths and apparent 
column densities, N a (Savage & Sembach 1992). These measurements are presented in Table 4. 
We also list our adopted central wavelengths and oscillator strengths (all from Morton 2003) for the 
transitions, the velocity range over which the interstellar absorption lines were integrated, and an 
empirical estimate of the signal to noise ratio per resolution element in the region about each line. 
The quoted measurement errors include contributions from continuum placement uncertainties and 
the effects of a 2% error in the flux zero level (c.f., Sembach & Savage 1992). 

Of the metal species considered in this work, only S II is detected in the STIS data. The 
absorption line profiles oftheS II triplet at 1250.584, 1253.811, and 1259.519 A from STIS E140M 
observations of vZ 1128 are shown in Figure 2; also shown is the profile of the S III 1012.495 A 
line from the FUSE observations of this star (Paper I). To calculate the limiting equivalent widths 
and apparent column densities of S I and Fe I, we assume the gas traced by these species has an 
intrinsic width of ~ 40 km s _1 (FWHM), which is the value derived from a single Gaussian profile 
fit to the S II 1250.584 A transition. All limits given in this work are 3a. 

The integrated apparent column densities of the weaker two transitions at 1250.584 and 
1253.811 A, which have /-values that are different by a factor of two, differ by ~ 0.01 dex, 
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with the weaker transition giving a higher N a . The strongest line, with an /-value three times that 
of the weakest, gives a value of iV a that is -0.05 dex lower than the weakest line. The progression 
of decreasing iV a with increasing / suggests these lines may contain unresolved saturated structure 
(Savage & Sembach 1991). 

Figure 3 shows a comparison of the apparent column density profiles as a function of velocity, 
N a (v), for the three transitions. In each of the panels, the weakest S II transition at 1250.584 A 
is shown as the thin black histogram, while the two stronger transitions are plotted as thick gray 
histograms. The N a (v) profiles in Figure 3 do not show the classical symptoms of unresolved 
saturated structure, with the highest peaks of the strong line profiles suppressed compared with the 
weak line profiles. The central dip near -ulsr ~ — 20 km s _1 in the S II 1259 A profile may be 
indicating the presence of a blended, but saturated, component at these velocities. 

Because of the possible presence of saturated structure in these profiles, we apply a correction 
to the data following the methodology outlined in Savage & Sembach (1991). We derive this 
correction separately for the line pairs A1250+A1253 and A1253+A1259 and average the results. 
Our adopted total column density is logiV(S II) = 15.28 ± 0.02. Independent support for this 
column density comes from a curve of growth analysis. A single component curve of growth fit 
to the equivalent widths of these three transitions gives logiV(S II) = 15.28lo!o2 f° r a Doppler 
parameter of b = 19.5 ± 1.2 km s _1 . 

Table 5 lists the final column densities adopted in this work, including those of S III, Fe II, 
Fe III, and O VI and of the limits to H 2 , S IV, and S VI from our analysis of the FUSE data in 
Paper I. 



3.2. Pulsar Dispersion Measurements and the Electron Column Density 

3.2.1. Radio Pulsar Dispersion Measurements 

The electron column density toward M 3 is derived from the dispersion measures determined 
by Ransom et al. (2004) for three binary millisecond pulsars in M 3. These pulsars were discovered 
with the Arecibo 305-m radio telescope at 20 cm using the Wideband Arecibo Pulsar Processors 
and the search algorithms described by Ransom, Cordes, & Eikenberry (2003) and Ransom, Eiken- 
berry, & Middleditch (2002). 

The average dispersion measure toward the three pulsars M 3 A, M 3B, and M 3D is DM = 
26.33 ± 0.15 pc cm -3 (standard deviation). We do not include the unconfirmed pulsar M 3C in 
this average, although it gives a consistent DM (Ransom et al. 2004). Typical uncertainties in 
the individual measurements are ~ 0.1 pc cm~ 3 . The average dispersion measure corresponds to 
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an electron column density of log N(e ) = 19.912 ± 0.002. This is the total electron column, 
including contributions from the WIM and the HIM along this sight line. 



3.2.2. Correcting for the Hot ISM 

Applying equation (3) to the 10 kpc sight line toward vZ 1 128 (Paper I), which is very close to 
the north Galactic pole, yields N(c~)him ~ 1.2 x 10 19 cm" 2 . This estimate assumes an interstellar 
gas-phase abundance A(0) = 4.08 x 10~ 4 (Andre et al. 2003), an ionization fraction x(0 +5 ) = 0.2 
(Sutherland & Dopita 1993; see Savage et al. 2003), and an O VI column JV(0 VI) = 3.1 x 10 14 
cm" 2 (Paper I). We adopt (n H ) c ~ 2 x 10~ 4 cm" 3 as the coronal density for the second term in 
brackets in equation (3). We have assumed a medium with fully-ionized hydrogen and helium, 
Vhim = 0.83. 

We will adopt a 50% uncertainty in calculations of the HIM column using the parameters 
for equation (3) adopted above. We have assumed x(0 +5 ) equal to its maximum in collisional 
ionization equilibrium models (Sutherland & Dopita 1993). While much of the O VI may in fact 
reside in regions where x(0 +5 ) 0.2, the O VI ionization fraction in the transition temperature 
gas along extended sight lines through the thick disk of the Galaxy is uncertain. We've assumed 
a single, constant density component for the 10 6 K material in an effort to approximate a Galactic 
corona. Yao & Wang (2005) have estimated the density distribution of this hot component of the 
ISM using X-ray absorption measurements of Ne IX, O VII and O VIII toward several Galactic 
targets, mostly in the Galactic Center region. They present two model density distributions based 
on their data: a plane-parallel, exponential disk distribution and a spherical distribution about the 
Galactic Center. These models give estimates of N(e~)niM ~ 2.9 x 10 19 cm" 2 and ~ 1.2 x 10 19 
cm" 2 , respectively, for the sight line to vZ 1 128 (with uncertainties > 30 to 50%). However, these 
results are likely to be strongly influenced by hot gas in the inner Galaxy. It is not clear that the 
Yao & Wang density distributions are appropriate for the hot ISM at the solar circle. We will adopt 
the predictions of equation (3) as calculated in the preceding paragraph for the HIM column and 
assume a 50% uncertainty. 



3.3. Total Abundances 

Table 6 summarizes the results of our analysis of the M 3 sight line, following the method 
described in §2. The total column density of sulfur toward vZ 1128 is logiV(S) = 15.34 ± 0.02; 
the contributions from S I, S IV and S VI are negligible [with combined 3a limits less than a few 
percent of iV(S)]. The total column of iron is log iV(Fe) = 14.95 ± 0.04, where we have assumed 
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the contributions from Fe IV and higher ionization states in WNM and WIM gas are negligible (as 
they are for ions of sulfur). 

In Table 6 we include the values of A(S) and A(Fe) for both the minimum and maximum 
helium ionization assumptions. We give the values adopted for the helium correction, r] WIM and 
the values of N(e~) associated with the WIM and the HIM. Also included in Table 6 are the 
column of hydrogen (neutral+ionized) and the fraction of hydrogen found in the ionized phase. 
Those values with a "Total" subscript refer to measurements that include the contribution from the 
HIM. Thus, the hydrogen column N(R) Tota i includes contributions from the WNM, WIM, and 
HIM, whereas iV(H) includes contributions from the WNM and WIM only. As stated above, the 
reason for this distinction is that we have not measured metals associated with the hot phase with 
the exception of O VI. 



4. Discussion and Conclusions 

We have derived the average gas-phase interstellar abundances in a way that minimizes ion- 
ization uncertainties along the sight line to the halo globular cluster M 3. The abundance of sulfur 
toward M 3 is log A(S) = —4.87 ± 0.03. Because sulfur is not believed to be significantly in- 
corporated into dust grains, A(S) represents the total average metal abundance of the neutral and 
ionized ISM along this sight line. The sight line to M 3 probes warm material associated with the 
thin interstellar disk and the extended (thick disk) distribution of neutral and ionized gas (Paper I). 

Typically, the solar system abundance for sulfur adopted by ISM studies is logA(S) = 
—4.73 (Savage & Sembach 1996), an average of the discrepant photospheric, log A(S) = — 4.67± 
0.11, and meteoritic, log A(S) Q = -4.80 ± 0.06, values from Grevesse & Sauval (1998). Thus, 
our total sulfur abundance along the sight line to M 3 is —0.14 ± 0.03 dex below, or ~ 3/4 of, 
the commonly-adopted solar value. 2 Recent revisions to the solar composition, however, are in 
very good agreement with our new measurement. Asplund, Grevesse, & Sauval (2005) derive 
log A(S) = —4.86 ± 0.05 for the solar photosphere, while incorporating the meteoritic data from 
Lodders (2003) with their results gives a meteoritic abundance of log A(S) = —4.84 ± 0.04. 3 It is 
somewhat disconcerting for ISM studies that the recommended solar abundances are in such a state 
of flux. However, the excellent agreement between our determination of A(S) toward vZ 1 128 and 



2 In the remaining discussion we adopt the results of the minimum helium ionization case, our preferred model. For 
simplicity we do not explicitly refer to the ~ 0.03 dex systematic uncertainty associated with this choice of model. 

3 Because H is deficient in meteorites, the metal abundances in meteorites are derived relative to Si. The absolute 
meteoritic abundances are then calculated from an adopted solar photospheric Si/H ratio. In this case, Asplund et al. 
(2005) use log A(Si) Q = 7.51 ± 0.04. 
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the most recent solar system estimates is encouraging. 

In contrast to sulfur, iron is heavily incorporated into grains. Our iron abundance measurement 
is log A(Fe) = —5.26 ± 0.04. Most of the difference between the gas-phase interstellar abundance 
and the solar system abundance, log A(Fe) Q = —4.55 ± 0.03 (Asplund et al. 2005), is due to the 
incorporation of Fe into dust in both the WNM and WIM. 

The ratio of sulfur to hydrogen derived in the usual way toward vZ 1 128 is log iV(S II) /N(H I) = 
—4.70 ± 0.04, consistent with the Grevesse & Sauval (1998) solar system value and +0.17 dex 
higher than the total abundance derived above using our method of correcting for ionized gas 
contributions. The neglect of ionization corrections causes a systematic error that is significantly 
higher than the estimated statistical uncertainties. Similarly, we derive logiV(Fe Il)/A r (H I) = 
—5.18 ± 0.06, implying a systematic error of +0.09 dex. The abundances derived by comparing 
singly-ionized sulfur and iron with H I both overestimate the true gas-phase abundances by signif- 
icant amounts. We note that although the direct comparison of two singly-ionized species reduces 
the impact of ionization uncertainties, they are still not negligible compared with typical statistical 
uncertainties. For the sight line toward M 3 we find log iV(Fe) /N(S) = —0.40+0.05, while simply 
comparing the column densities of singly-ionized atoms gives logiV(Fe Il)/iV(S II) = —0.48 ± 
0.05. For comparison, the solar system ratio from Asplund et al. (2005) is logiV(Fe)/iV(S) = 
+0.31, while the standardly-adopted value in ISM studies has been +0.23 (see Savage & Sembach 
1996). 

The overestimate of gas-phase abundances derived by comparing singly-ionized species with 
neutral hydrogen is not unique to this sight line. While the fraction iV(H ll)Totai/N(R) Tota i ps 
0.45 ± 0.04 for the vZ 1128 sight line is the highest yet reported (Reynolds 1993), the WIM will 
contribute to the singly-ionized metal column densities along all high-latitude sight lines. Several 
groups have derived near-solar abundances of sulfur comparing S II with H I along high-latitude 
sight lines with H I column densities similar to that toward vZ 1 128 in M 3 (e.g., Howk et al. 1999; 
Spitzer & Fitzpatrick 1993). It is likely that the abundances derived in these studies suffer from 
the systematic effects of ionized gas as well. 

It may be possible to model the contribution from the ionized gas along a sight line and correct 
for its systematic effects on abundance studies. Using the WIM models constructed by Domgorgen 
& Mathis (1994), Sembach & Savage (1996) estimated that the magnitude of ionization uncertain- 
ties for derived gas-phase elemental abundances should be less than ~ 0.1 to 0.2 dex along most 
sight lines through the Galactic halo, consistent with the values derived in this work. 

Unless A(S) varies significantly within the Galaxy, however, correcting for ionization effects 
may not be so straightforward. Howk et al. (1999), in their detailed study of the sight line to 
fx Columbae, derived log iV(S Il)/iV(H I) = —4.65 ± 0.02; they used photoionization models to 
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argue that the sulfur abundance along this sight line was near their adopted solar system value of 
log A(S) = —4.73 (as recommended by Savage & Sembach 1996). In this case, the star is much 
closer to the Sun than vZ 1 128, and Howk et al. assumed an H II region about the star dominated 
the ionized gas content of the sight line. Even adopting the maximum ionization correction derived 
by Howk et al. (1999), however, gives a sulfur abundance ps 0.2 dex higher than that derived toward 
vZ 1 128 in this work. 

Sembach et al. (2000) produced detailed models for the ionization of the WIM by OB stars 
and used the models to estimate the ionization corrections required for the sight line to the distant 
halo star HD 93521 (Spitzer & Fitzpatrick 1993). Their models predict an ionization correction 
of ps —0.2 dex to A(S) derived by Spitzer & Fitzpatrick (1993), giving A(S) consistent with that 
derived in this work. However, applying their method to the sight line toward vZ 1128 yields an 
abundance log A(S) = — 5.3, or ps —0.4 dex from the value derived in this work. 

These examples call into question our ability to effectively model the effects of WIM gas 
more generally. Indeed, models designed to match emission line diagnostics may not be good 
descriptors of the absorption line diagnostics of the WIM, which are more sensitive to low-density 
gas. 

We can state with confidence that abundance studies of almost all high-latitude, low hydrogen 
column density [log iV(H I) < 20] sight lines are affected by the presence of ionized gas when 
comparing singly-ionized metal species with H I. The degree to which an individual sight line is 
contaminated likely depends on the fraction of hydrogen along the sight line that is ionized. The 
work of Reynolds (1993) suggests that this fraction can range from ~ 20% to ~ 50%. We note that 
some species can be reliably compared if they have similar ionization properties; for example, the 
ratio O I/H I should not be affected by ionization effects like those discussed here given the strong 
charge exchange reaction that locks the ionization state of these species together. Unfortunately, 
the far-UV O I lines toward vZ 1 128 are strongly saturated while the weak O I 1355 A transition 
is not detected. 

We will apply the technique outlined here to more sight lines and to an extended array of 
metal species in the future to study this source of systematic uncertainties in regions of varying 
physical conditions and fractional ionization. At least three other globular clusters contain both 
radio pulsars and stars with sufficient UV flux to allow measurements like those presented in this 
work. 

We thank S. Ransom for discussions regarding the dispersion measures toward pulsars in 
M 3. Support for proposal numbers HST-GO-9150.03-A and HST-GO-9410.03-A was provided 
by NASA through grants from the Space Telescope Science Institute, which is operated by the 
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Table 1 . Values of Helium Correction Factor rj 



x(R+) x(Re + )/x(R+) x(He +2 ) 77 



Minimum Helium Ionization 21 



0.67 


0.15 


0.00 


0.99 


0.80 


0.15 


0.00 


0.99 


1.00 


0.15 


0.00 


0.99 


0.67 


0.27 


0.00 


0.97 


0.80 


0.27 


0.00 


0.97 


1.00 


0.27 


0.00 


0.97 


Maximum Helium Ionization 13 


0.67 


0.00 


1.00 


0.77 


0.80 


0.00 


1.00 


0.80 


1.00 


0.00 


1.00 


0.83 c 


0.67 


0.15 


0.78 


0.80 


0.80 


0.15 


0.81 


0.82 


1.00 


0.15 


0.85 


0.84 


0.67 


0.27 


0.60 


0.83 


0.80 


0.27 


0.66 


0.84 


1.00 


0.27 


0.73 


0.85 



a The minimum helium ionization case as- 
sumes x(He +2 ) = 0. 

b The maximum helium ionization case as- 
sumes x(Re +2 ) = 1 - x(Ke + ). 

c This case is appropriate for HIM gas. 
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Table 2. STIS Observations of vZ 1 128 a 



Archive ID 


Grating 


t e xp [ksec] 


R = A/AA 


FWHM [km/s] 


A Range [A] 


O6F502010 


E140M 


2.0 


45000 


6.5 


1140-1735 


O6F502020 


E140M 


10.7 


45000 


6.5 


1140-1735 


O6F501010 


E230M 


2.0 


30000 


10 


1575-2380 


O6F501020 


E230M 


8.0 


30000 


10 


1575-2380 


O6F501030 


E230M 


2.7 


30000 


10 


2305-3110 



a All observations were taken through the 0/2 x 0'/06 apertures on 2002 August 
16-17. 



Table 3. HI Observations in the Direction of vZ 1 128 



Source 


Beam 


logiV(H I) a 


Source b 


NRAO 21 cm 


21' 


20.058 ±0.003 


1 


LDS 21 cm 


35' 


20.044 ±0.005 


1 


LDS 21 cm 


35' 


20.017 ±0.006 


2 


STIS Lyman- a 




19.98 ±0.03 


3 



a The uncertainties quoted for the 21 cm observa- 
tions do not contain a contribution from errors in 
the stray radiation or baseline corrections. These are 
likely of order 0.04 dex for this sight line. 

b Sources: (1) This work, through direct integration 
of the brightness temperature profile; (2) Wakker et 
al. 2003 through a Gaussian decomposition of the 
profile; (3) This work, through a fit to the Lyman-a 
profile (see text). 
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Table 4. Selected Interstellar Absorption Line Measurements toward vZ 1 128 



Species 


A c [A] b 


logA/ c 


v _, v + d 


W x [mA] 


logiV a 


S/N 


S I 


1425.030 


2.437 




< 15 


< 12.7 


16 


S II 


1250.584 


0.834 


-50,35 


102 ±4 


15.25 ±0.02 


19 


S II 


1253.811 


1.135 


-75,35 


162 ±5 


15.24 ±0.02 


19 


S II 


1259.519 


1.311 


-75,35 


190 ±6 


15.20 ±0.03 


15 


Fe I 


2167.453 


2.512 




< 21 


< 12.5 


18 



Measurements from STIS E140M and E230M observations of vZ 1 128. 
All limits are 3a and uncertainties include contributions from photon statis- 
tics and an assumed zero-point uncertainty of 2% of the local continuum. 
Signal-to-noise ratios, S/N, determined empirically following Sembach & 
Savage 1992. 

b Rest wavelength of the transition, from the compilation of Morton 2003. 

c Oscillator strength, /, given as log Xf, from the compilation of Morton 
2003. 

d LSR velocity range of integration for quoted values of equivalent width 
and apparent column density. 
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Table 5. Adopted Column Densities Toward M 3 



Species 


logiV 


Ref. a 


HI 


19.98 ±0.03 


1 


e~ 


19.912 ±0.002 


2 


H 2 


< 14.35 (3a) 


1 


VI 


14.49 ± 0.03 


1 


S I 


< 12.7 (3d) 


3 


S II 


15.28 ±0.02 


3 


S III 


14.47 ±0.03 


1 


S IV 


< 13.7(3(7) 


1 


S VI 


< 13.4(3(7) 


1 


Fe I 


< 12.5 (3a) 


3 


Fe II 


14.80 ±0.05 


1 


Fe III 


14.42 ± 0.05 


1 



a References: (1) Howk et al. 
2003 (Paper I); (2) Ransom et al. 
2004; (3) this work. 
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Table 6. Derived Interstellar Parameters Toward M 3 



Helium Ionization 



On anti fv 
u chili lv 


lVTi n i tn 1 1 tn a 

1VA1111111 LI 1 1 1 


Maximum* 3 


VWIM ° 


0.98 ±0.01 


u.oi in u 


nzL 


log Me - ) 


19.912 ± 


0.002 




\ogN(e') HIM 


19.10 ± 


0.12 




\ogN(e~) WIM 


19.84 ± 


0.03 




logiV(H R)him 


19.01 ± 


0.18 




logiV(H \\) WIM 


19.83 ±0.04 


19.75 ±0. 


.04 


logiV(H) d 


20.22 ±0.02 


20.19 ±0. 


03 


\ogN(Yi) Tota f 


20.24 ±0.03 


20.21 ±0. 


03 


N(H Il)H, /M /iV(H) d 


0.42 ±0.04 


0.37 ±0. 


.04 


N(U ll)Total/N(R) Total C 


0.45 ±0.05 


0.41 ±0. 


.06 


logW(S) 


15.34 ± 


0.02 




logiV(Fe) 


14.95 ± 


0.04 




logA(S) 


-4.87 ±0.03 - 


-4.84 ±0. 


03 


logA(Fe) 


-5.27 ±0.05 


-5.23 ±0. 


.05 



a The minimum helium ionization case assumes 
x(iie +2 ) WI M = 0. We prefer this assumption over the 
maximum helium ionization case. 

b The maximum helium ionization case assumes 
x(Re +2 ) WIM = 1 - x(He + ) WIM . 

c The values of r] WIM are derived assuming 0.67 < 
x(R + )wim < 1.0 and x(Re + ) WIM < 0.27 x(R+) WIM . 
The uncertainties quoted encompass the range of values de- 
rived using these constraints and the assumptions regarding 
x(Re +2 ) WIM appropriate for each model. 

d These quantities refer only to warm gas, including the 
WNM and WIM. They exclude ionized gas associated with 
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the HIM. 

c The quantities with the subscript "Total" include con- 
tributions from both warm and hot gas, i.e., they include 
ionized gas associated with the HIM. 
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Fig. 1. — Two views of the STIS E140M spectrum of the Lyman-a; absorption profile toward 
vZ 1128. Both stellar and interstellar absorption contribute to this profile. The estimated stellar 
profile is shown as the thin blue line and has been shifted to vlsr = —145 km s _1 to match the 
observed positions of stellar absorption lines in the STIS spectrum. The best fit to the interstel- 
lar and stellar absorption profile is shown in red, corresponding to an interstellar column density 
log iV(H I) = 19.98 ± 0.03. The green line shows the best fit profile adopting the H I column den- 
sity derived from Green Bank 140-ft telescope observations of 21-cm emission in this direction, 
logiV(H I) = 20.04 ± 0.01 (see Paper I). The sharp emission line in the center of the Lyman-a 
absorption trough is geocoronal emission. The strong stellar line near 1206 A is stellar Si III. 
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Fig. 2.— Absorption line profiles of S II A1250.584, 1253.811, and 1259.519 from STIS E140M 
observations and S III A1012.495 from FUSE (see Paper I). We suggested in Paper I that the 
material at v lsr ~ —30 km s _1 was likely associated with the Galactic thick disk given the higher 
degree of ionization compared with the material at vlsr ~ —4 km s _1 . The STIS data have a 
resolution of ~ 6.5 km s _1 , while the FUSE data have a resolution ~ 20 km s _1 . 
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Fig. 3. — Comparison of the apparent column density profiles for S II A1250.584, 1253.811, and 
1259.519 from STIS E140M observations. The weaker two lines (at 1250.584 and 1253.811 A) 
are compared in the top panel; the strongest line (1259.519 A) is compared with the weakest in the 
lower panel. The lack of agreement between these profiles indicates the stronger line contains a 
small amount of unresolved saturated structure. 



